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Abstract. 

The Bogoliubov model for weakly interacting Bose gas is extended to Bose-Einstein 
condensation (BEC) of spin-S atoms in a magnetic field. Equation for the vectorial 
order parameter valid at temperature T — > is derived and its particular solution is 
found. This solution corresponds to the formation of BEC of atoms with a definite spin 
projection onto direction of a magnetic field. We study the thermodynamic stability 
of the found solution and obtain the expressions for low-lying collective modes. 
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1. Introduction 

After the first remarkable experiments concerning the observation of BEC in dilute gases 
of alkali atoms such as ®^Rb [1], ^^Na [2], and ^Li [3] the interest to this phenomenon has 
revived [H [5] . Later on, BEC has been also obtained in other atomic species: atomic 
hydrogen [6], metastable ^He [7], and ^^K |8]. The experimental realization of BEC has 
become possible due to the progress of laser cooling and trapping techniques [9]. The 
carried out experiments have proved many predictions of the microscopic theory for 
weakly interacting Bose gas, which originates from the pioneering work of Bogoliubov 
[lOj . Bogoliubov's theory has become almost the first theory in which it was necessary to 
move essentially from the methods of standard perturbative approach while describing 
the interaction effects. However, this theory, in its original formulation, did not take 
into account the internal degrees of freedom of atoms. The effect of spin degrees of 
freedom for weakly interacting Bose gas (spinor BEC) has been studied in |llj-[19j. 

The realization of optical trapping for atomic condensate [20] has stimulated 
theoretical interest to spinor BEC. Bose condensation in a weakly interacting gas of 
bosonic atoms has been studied theoretically by many authors both for spin-1 p!2]-[T7] 
and spin-2 [18], [19] bosons. These investigations are based on the effective interaction 
Hamiltonians of two bosons, in which the interaction is characterized by a definite 
number of interaction constants - s-wave scattering lengths. The number of scattering 
lengths is determined by the total spin of two interacting bosons taking into account the 
symmetry properties of their wave function. For example, in case of spin-1 atoms the 
interaction Hamiltonian contains two interaction constants [12]- [17], in case of spin-2 
atoms there are three interaction constants [18], [19]. Thus, as the spin value of atoms 
grows, the number of constants, which characterize the interaction of two bosons, is 
increased. Note that in the mentioned effective Hamiltonians it is difficult to interpret 
the physical nature of the separated term of non-relativistic interaction not associated 
with neither potential nor spin-exchange interactions (see e.g. [18]). 

In this paper we study a weakly interacting Bose gas of particles with arbitrary 
integer spin in a magnetic field (see also [H]). We start from the microscopic 
interaction Hamiltonian for two spin-S* bosons. This Hamiltonian is specified by two 
functions, which describe potential and spin-exchange interactions of spin-S* atoms. 
According to general rules of quantum mechanics, we pass from the pairwise interaction 
of two bosons to the standard expression for binary interaction of arbitrary number of 
bosons in the second quantization representation. By solving the multichannel scattering 
problem for the considered Hamiltonian we could find, in principle, all scattering lengths 
in terms of the functions characterizing the potential and spin-exchange interactions. 
Thereby, it would be possible to obtain the Hamiltonians analogous to the above 
mentioned effective interaction Hamiltonians (see e.g. [E]). However, the use of the 
microscopic Hamiltonian enables to restrict ourself by two interaction constants even in 
the case of arbitrary spin when studying the ground state, stability, and excitations in 
a weakly interacting gas in the presence of BEC. 
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2. Method of quasiaverages and the model with a separated condensate 

To describe the system with a spontaneously broken symmetry we address to the method 
of quasiaverages [21], [22]. According to this method the Gibbs statistical operator is 
modified so that it possesses the symmetry of degenerate state. This modification is 
usually done by introducing the infinitesimal "source" uF {y 0) into the Gibbs 
exponent, which has the symmetry of phase under consideration. Then, the average 
value of any physical quantity A is defined as 

-< A >~= lim lim Ttw^A, (1) 

where the Gibbs statistical operator Wy has the form 

Wy = exp (n^ - (3{H - /iiV + uF)) . (2) 

Here (3 = 1/T, /i are the reciprocal temperature and chemical potential respectively 
and H, N are the system Hamiltonian and the particle number operator. The 
thermodynamic potential fli, being a function of thermodynamic parameters (3, fi is 
found from the normalization condition Trt&j, = 1. Notice that the limits in ([1]) are not 
permutable. 

Consider a gas of condensed bosonic atoms with spin S. The formation of a 
condensate is accompanied by the gauge symmetry breaking and, therefore, in order 
to remove this kind of degeneracy we should choose the "source" uF in ([2]) such that 
[wjy, A^] 7^ (A^ is the generator of phase transformation), 

vF = Vo.j c?'x(V^i(x) + V;,(x)), (3) 

where '?/'^(x), '?/'a(x) are the creation and annihilation operators with index a taking 
25* + 1 values (the summation over repeated indices is assumed). Then, according to 
([I]), (El), -< ■?/'a(x) V~^l'^ -< I that corresponds to the formation of atomic 

condensate with momenta p = 0. The order parameter "ifa = V'^l'^ -< a^a >- is called 
the condensate wave function. 

The method of quasiaverages and the spatial correlation decay principle enable to 
justify the replacement of creation and annihilation operators of atoms with momentum 
p = by c-numbers, doa, Oq^ — » VV'^a, VV'^a [2l]-[23] (the condensate separation 
procedure). 

The basic statement of the method of quasiaverages applied to the description of 
BEG consists in the following [21]- [23]: the Gibbs statistical operator is replaced by 

w{^) = exp (fi(^) - /?(^(^) - /iiV(^))), (4) 

where \Ef = {\E'q,, \E'* } is found from the following equation: 

dm) n (^^ 
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3. The ground state of spin-S* condensate in a magnetic field 

In this section we study one of the possible ground states of spin-S BEC in a magnetic 
field. In doing so, we start from the Hamiltonian Ti = H — /iiV, which determines the 
Gibbs statistical operator (jlj) and has the following form: 

Ti. = Ti-o + Ti-p + Ti-e, (6) 



where 



"^P = X! ^(Pl3)<^Pi+P2,P3+P4^pia^P2/3^P3a^P4/35 (8) 

P1...P4 

~ ^Tf7 X/ '^(Pl3)'^Pi+P2,P3+P4^pia^P2/3^f*7^/3p^P37^P4P- (9) 
P1...P4 

Here Sq/j are the spin matrices, t/(pi3), ^(pis) (pis = Pi — Ps) are the Fourier 
transforms of the amplitudes of potential and spin-exchange interactions respectively, 
and h = gtl/S {g is the Bohr magneton, and H is an external magnetic field). 
For our next calculations it is convenient to introduce the so-called ladder operators 
S± = Sx iSy. Then, their nonzero matrix elements in the representation, where 5*2 is 
a diagonal matrix, {a\Sz\a') = aSaa', have the form 



a + l\S+\a) = JS{S + 1) -a{a + l), (10) 



{a - l\S^\a) = ^S{S +1) -a{a-l). 

Now we separate the p = components Soa in the Hamiltonian (the replacement 
of O'Oa by c -numbers, doa VV'^a) and keep the terms only up to second order in 
dpa- We omit the higher order terms, since they should be taken into account only 
when examining the interaction between quasiparticles, which we will introduce in the 
next section. As a result the Hamiltonian takes the form Ti. ~ T-C^^^ + Ti^"^^ . The explicit 
expression for 7^^°\ which contains only c -numbers \E'q, reads 



V 

where 



(^*^)^ + ^^(^*S^)^-h^*S^-/i^*^, (11) 



= ^*S^ = ^ISo^p^p (12) 

The explicit form for Ti*^^^ will be written in the next section. 

Next, making use the normalization condition Trw) = 1, we find immediately the 
thermodynamic potential density uj = QT/V in the leading approximation (neglect of 
quasiparticles; T — > 0) of the model for weakly interacting Bose gas, 

UJ = ^^(^*^)2 + :^(^*s^)2 _ h^*s^ - /i^*^. (13) 

Therefore, Eq. ([5]) for \1'q, takes the form 

- f/(0)(^*^)^, - J(0)(^*S^)S,^^^ + hS,^^^ = 0. 
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If to introduce the normahzed spin functions (a, = a/^Cc^ where n = is the 
condensate density and CaCa — 1' then the latter equation is written as 

- nUm^ - nJ{0){CSC)Sai3Ci3 + hSapCf, = 0- (14) 

Assuming the vector h directed along ^-axis (h = {0,0, h)), its solution Ci™^ being an 
eigenfunction of S^, iSz)ai3Cjr^ = ^C^\ has the form 

= Sam. (15) 

Next, taking into account that Sz is a diagonal matrix, whereas S± have no diagonal 

matrix elements in the considered representation of spin matrices, one finds from Eq. 



f/(0)+m2J(0)- ^ ' 

The obtained formulae (ITSll . (fT6l) result in the following expression for the 



thermodynamic potential density: 
_ 1 (/i + mhf 

2f/(0)+mV(0)' ^ ' 

We are now in a position to study the stability of possible ground states (fT5l) . In the 
considered approximation, the thermodynamic potential of the normal state is zero (the 
order parameter '^^ vanishes). Therefore, for the stability of the ground state under 
consideration, the density of thermodynamic potential must be negative, uj < and, 
consequently, according to ( |T71) . we can write the necessary condition of thermodynamic 
stability, 

f/(0) + mV(0) > 0. (18) 

Let us find now such spin projections m, which correspond to the minimum of potential 
f|T7|) . For simplicity, we study the case oi h = (or sufficiently weak h). Then, 

1 ^„ 
2[/(0)+mV(0) 

As it can be easily seen that in contrast to usual Bogoliubov's theory, in which U{0) > 
(the necessary condition of stability), the negative values of U{0) are also permissible. 
Therefore, we have the following three situations: 

1) U{0) > 0, J(0) > 0. In this case the requirement (ITSl) is automatically satisfied. 
The density of thermodynamic potential (ITTll has a minimum at m = in which 
LJ = — /x^/2[/(0). We call this case as antiferromagnetic ordering. 

2) U{0) > 0, J(0) < but such that the requirement ( fTSl) should be satisfied. The 
minimum of u is reached for mmin = ^{fnc — 1), where 

mc= ("tII) ^ ' m,<S + l (19) 

(the square brackets denote an integer part). This case corresponds to ferromagnetic 
ordering. 
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3) U{0) < 0, J(0) > but again, such that U{0) +m'^J{0) > 0. Here the minimum 
of u is given by the spin projections rrimin = ±(^c + 1)? where nic is also defined by f|T9l) 
but with f/(0) < 0, J(0) > 0. This case also corresponds to ferromagnetic ordering. 



4. Low-lying collective modes 



In this section we obtain the excitation spectra of spin-S* BEC by employing the well- 
known diagonalization procedure (Bogoliubov's u — v transformations [10]) for the 
Hamiltonian quadratic in creation and annihilation operators. Note that the excitation 
spectra can also be found as a result of the linearization of the Gross-Pitaevskii equation 
[211 [25] for the condensate wave function (see e.g. |12]). 

The part of the spin-exchange interaction Hamiltonian ([9]), which is quadratic in 



"pa ; 



(p 7^ 0), has the form 



iif) = J(0)vI/*Svl>^atSap+-^ J(p) [(aJ.Sv|/)(alpSvI/) + (at,Sv|/)(vl>*Sap) + h.c] , (20) 
p ^ p 

where we have used the notations (fT2l) . Taking into account that Sx = 1(5'+ + S*-), 
Sy = —^{S+ — S-) and bearing in mind ffTOl) for non-zero matrix elements of Sj- as well 
as the explicit form of the condensate wave function \Ef2"'' = V^^am, one gets 



(vi/*SM/)(atSap) 
(at,SvI/)(aLpSvl/) 



nm > aa'a 



po"'pa) 



p»7i— 1 j ■> 



n 



(«pS^)(^*Sap) = nm a^^ap^ + - (S'_„a^„„iap™_i + S^al^_^_^apm+l 
where the following notation has been introduced: 

Srn. ~ 



S{S +1) -m{m + l). 
Hence, Ti^^-* takes the form 

H^^^ = J{0)nmY^ [(m - l)aj,„_iapm_i + mal^aprn + {m + l)a^„+iapm+i 



+J{0)nmJ2 



po^-pa ■ r) 
p,a ^ p 



m 



^pm^— pm ~^ ^^pm^pm ~^ pm^pm 



(^pm-l^-pm+l ~l~ ^-pm+l^^pm— ly 



+|E^(p) 
^ p 



r(2 -t " _|_ Q2 - 



2 ;it A 

m^pm+l^pm+l 



a 7^ m — 1, m, m + 1. 



In this formula, the summation index a in the second term takes all values of spin 
projections except m — 1, m, and m + 1 (these three projections we have separated off 
and written them as the first term in Ti.^^)- In a similar manner, keeping the terms only 
of second order in apa, one finds according to ([7]), ([8]) 



E/ (^pm-l^pm-l + ^pm^pm + ^pm+l^pm+ly 
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-h I (m - l)a|,„_iapm_i + maj^^ap^ + (m + l)aj, 



+ ^ (^p — /i — ah) dp^cipa, a 7^ m — 1, m, m + 1, 



and 



^pm+l^pm+l 



^pQ^pa 



(21) 



(22) 



+ „ E] ^(P) ^pm^-pm + 2ap^apm, + 
^ P 



a m — 1, m, m + 1. 



When obtaining fl2T]) we have employed the fact that h is directed along 2;-axis, 
h = {0,0, h). Next, using Eq. ( fT6l) to eliminate the chemical potential /z in ( 12T1) . 
we recast the total Hamiltonian Tt!^^) = + Ti^^^ + that is quadratic in creation 
and annihilation operators in the following form: 

where 



(23) 



^ = E [^p — J{0)nm{m — a) + h{m — a)] 0)^^(1^^, a ^ m — 1, m, m + 1, (24) 



p,a 



1 



^pm^— pm ~^ ^pm,^— pm 



(25) 



n^'^\m - l,m + 1) = ^ [^p - /l + Pmip)] apm+lOpm+l 

p 

^pm— l^pm— 1 

P 

+ 51 '^'"(P) ^pm-l^-pm+1 + ^Pm-l^-pm+1 • (26) 
P 

The introduced quantities am{p), Pm{p), and gm{p) are given by 

Tl 

«m(p) = -J{p)SmS-m, (27) 

/3^(p) = -J(p)S'^ + nJ(0)m, (28) 
gUp)=n{U{p)+m'j{p)). (29) 

Now we are in a position to carry out the diagonalization procedure of the total 
Hamiltonian fl23|) quadratic in creation and annihilation operators. In this connection we 
note that the " Hamiltonians" (l2^ - fl26l) contain the creation and annihilation operators 
with not overlapping sets of indices a, m, m — 1, m + 1 {a ^ m — 1, m, m + 1). 
Therefore, we can perform their diagonalization independently. The evidence of this 
statement also follows from the fact that (1231) can be considered as the Hamiltonian of 
the system consisting of four kinds (m, m ± 1, a) of noninteracting particles. 

The "Hamiltonian" Ti'^^ has already a diagonal form with the following spectrum: 

^m,a{p) = £p — J{0)nm{m — a) + h{m — a). (30) 
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To carry out the diagonalization of 'H^'^\m — l,m + 1), we introduce the creation 
and annihilation operators bp^n+a (o" = il)) 

^pm+o- — '^m,o-(P)^pm+o- + '^m,o-(p)^-pm-o-; (31) 

in terms of which it has the diagonal form, 

n^^\m - l,m + 1) = J2^"^AP)Km+abpm+a + ^0, (32) 

p.cr 

where LUm,u{p) and Eq are the excitation spectrum and the ground state energy 
respectively. In order that the introduced operators ftpm+o-j ipm+a nieet the canonical 
commutation relations, the functions Um,aip), Vm,a{p) niust obey the relationships 

|Mm,a(p)P - \Vm,a{P)\'^ = 1' 
Um,a{p)Vm-a{-p) - Vm,a{p)Um,^a{-p) 

Next, noting that 



0. 



(33) 



[7i;(2)(m-l,m + l),, 



■'pm+tT 



and expressing the right-hand side of this formula through fepm+o-, one gets 



[?Tf(2)(m - l,m + l),apm+<^ 



^am(p)t;;; _^(-p) + 7m,aMm,a(p) j &pm+a 

- (a^(p)M;; _^(-p) +7„,<,w^,^(p)) 6lp„_^, 

where 7^,0- is defined by 

7m,a = £p + Pma{p) " f^/^, 

moreover Pma{p) depends on the product ma. On the other hand, the straightforward 
use of (pri) results in 



(34) 
(35) 



'lHP\m - l,m + l),aprri+<x = ^^m,f7(p)t^m -a(-p)&-pm-a - ^im,<7 (p)t^m,a (p) • (36) 

The comparison of flMl) with fl36|) gives the coupled equations for Um,a{p) and _^(— p), 

(7m,a - Wm,a(p)) Mm,a(p) + ^^(p)^^;^,-^ (-p) = 0, 

am{p)UmAP) + (7m -a + ^^m,a{P)) ^m,-a(-p) = 0- (37) 

The condition for the existence of non-trivial solutions to the coupled Eqs. fl37|) along 

with the definitions fl7n) - fl2I?l) . ([35D result in the following expression for the excitation 
spectrum: 



^m,a(p) = nma ( J(0) - ^J{p) 



± 



El + e^nJ{jp){S{S + l] 



/nJ(p)m 



1/2 



ah. 



{31 



At small p and J(0) < the obtained spectrum is real (the state is stable) if the spin 
projections m meet the inequality > m^, where 
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and the square brackets, as in (fT^ . are introduced to denote an integer part. 
The functions Um,a{p), ^^m,o-(p) are found from fl33|) . fl37j) and have the form 



UmAP) 



VmAP) 



m.a 



In fact, the functions ^^ ^-(p), fm,o-(p) do not depend on cr because, as it can be easily 
shown, the quantity tUm^o-lp) — 7m,cr is independent of a. The sign plus before the 
square root in ( l38l) corresponds (for a = 1) to the wave, which propagates in one 
direction, whereas the sign minus corresponds (for a = —1) to the wave propagating in 
opposite direction. Notice that the obtained spectrum, as well as (l30l) . contains only the 
spin-exchange interaction amplitude and does not depend on the potential interaction 
amplitude. 

When m = (the antiferromagnetic ordering), the excitation spectrum ( l38l) takes 
the form 



uj{p) = + epnJ{p)S{S +l)±h. 
In this case, for h = and p ^ we have 



Ap) = cp, c=^^mS{S+l). 
In ferromagnetic case (when m = S) the excitation spectrum is of the form 

uj{p) =ep+ '"'^^^^'^ (1 -a)+ nJ{0)Sa - ah. 

The similar mathematical manipulations with W'^^m) lead to another mode of 
excitation spectrum, which depends both on potential and spin-exchange interaction 
amplitudes. 



Umip) = + 2Epn{Uip) + m2 J(p)). (39) 

The stability region for this spectrum at small p is given by flTI?]) . 
The functions uAp) and fm(p) can be found immediately, 

«m(p) = , ^ , vAp) = I , ■ 

I^EpUJAP) I^EpUjAP) 

When J(p) = the excitation spectrum fl39|) coincides with the spectrum found by 
Bogoliubov |T0]. At small p, the spectrum has the following phonon behavior: 

o;^(p) = cp, c=^— (t/(0)+mV(0)). 

In this formula, as well as in (1391) . we have chosen the arithmetic value of the square 
root. 

In conclusion, we have studied BEC of atoms with arbitrary spin in a magnetic field 
on the basis of the model for weakly interacting Bose gas. We have derived the equation 
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which describes the ground state of spin-S* BEC at temperature T — and found its 
particular solution. This solution corresponds to the formation of BEC of spin-S* atoms 
with a definite spin projection m onto direction of a magnetic field that is also true for 
an ideal Bose gas [11]. The explicit expression for thermodynamic potential being a 
function of chemical potential and spin projection has been obtained. It generalizes the 
thermodynamic potential for weakly interacting Bose gas to the case when both potential 
and spin-exchange interactions act between bosons. The thermodynamic stability of 
the state under consideration has been studied and the spin projections, which give a 
minimum of thermodynamic potential, have been found. These projections are given 
by the integer part of the ratio of potential to spin-exchange interaction amplitude. 
The expressions for low-lying collective modes related to the solution f|T5|) have been 
obtained. Notice that Eq. ( fT4l) for order parameter has also other solutions different 
from ( ITSi) . Our present research deals with seeking such solutions. 
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